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1 Introduction 

Assume C R 3 is a C 2 - smooth bounded domain and Qt = H x (0,T). In 
this paper we investigate the boundary regularity of solutions to the principal 
system of magnetohydrodynamics (the MHD equations): 



d t v + (v ■ V)w - Av + Vp = rot H x H 
div v = 



in Q T , (1) 



d t H + rot rot H = rot(u x H) 
div H = 



in Q T . (2) 



Here unknowns are the velocity field v : Qt —> R 3 , pressure p : Qt ~ > R, and the 
magnetic field H : Qt — > K 3 - We impose on v and H the boundary conditions: 

v \anx(o,T) — 0, ^i?|anx(o,T) = 0, (rot -ff) T |a^ x ( ,t) = 0, (3) 

Here by v we denote the outer normal to dfl and H v = H ■ v, (rotiJ) r = 
rot-ff — v(rotH ■ v). We introduce the following definition: 

Definition: Assume T C <9fi. The functions (v,H,p) are called a boundary 
suitable weak solution to the system (1), (2) near IV = Tx (0, T) if 

1) v e £ 2 ,oo(Qt) n W^\Q T ) n wfi (Qt), 

8 ' 2 

^ei2,oo(QT)nw 2 1 ' (Q T ), 

2) P 6L3(g r )n^' 3 (Q r ), 

2 8 ' 2 

3) divt> = 0, div H = a.e. in Qt, 

4) = 0, H v \ga = in the sense of traces, 
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1 1 ^ 



J v(x, t) ■ w(x) 



dx 



5) for any w e £ 2 (£1) the functions 

and t n> J H(x,t) ■ w(x) dx 
n Q 
are continuous, 

6) (v, H) satisfy the following integral identities: for any t € [0, T] 

J v(x,t) ■ T](x,t) dx — J vq(x) ■ r)(x, 0) dx + 



+ J J ( - v ■ d t r)+ (Vv - v ® v + H (g> H) : Vi] - (p + ±\H \ 2 ) div rj) dxdt 
o n v ' 

for all r\ e Wl^iQt) such that r?|anx(o,t) = 0, 

/ H(x, t) ■ ip{x, t) dx - J H Q (x) ■ i>(x,0) dx + 
n n 

+ / / ( - H ■ d t tp + rot H ■ rot V - (v x H) ■ rot ip) dxdt = 0, 



o. 



for all V e ^^(Qt) such that 



aax(o,t) 



7) For every z = (a; ,io) € T T such that Q R (x ) x (i - R 2 ,ta) C Qt and 

= the following 



for any C G C^{B R {x ) x (t„ - R 2 , to]) such that §§ 
"local energy inequality near Tt" holds: 



sup / C(\v\ 2 + \H\ 2 ) dx + 
te(t a -R?,t a )n R (x ) v 7 

+ 2 / / c(|Vu| 2 + rot i^| 2 ) dxdt < 

t -R 2 n R (x ) 

< S I (\v\ 2 + \H\A(d t ( + A() dxdt + 

to , N (4) 

+ / / I |w| 2 +2p)v ■ V( dxdt + 
t -R 2 n R (x ) v ' 

to 

- 2 / / (H ® H) : V 2 C dxdt + 
t --R 2 n H (x ) 

to 

+ 2 / J* (v x i?)(VC x iJ) dxdi 
We remark also that the following identity holds 

to 

/ / (v x i?)(VC x H) dxdt = 

t -R 2 n R (x ) 

= J J (vVQ\H\ 2 dxdt - J J (v H)(H • VC) dxdt 

to-R 2 n R (x ) t -R 2 n R (x ) 



2 



Here L s j(Qt) is the anisotropic Lebesgue space equipped with the norm 

ll/lk,« T ) == (f a l/0M)| s dx) l/S dtf'\ 

and we use the following notation for the functional spaces: 

Wlf(Q T ) = Li(0,T;W}(O)) = {ue L,,i(Q T ) : V« G L sJ (Q T ) }, 
W 2 s f(Q T ) = {ue W]f(Q T ) : V 2 u, d t u e £ s j(Qt) }, 
^(S]) = {«e^(fi): w|an=0}, 
and the following notation for the norms: 

\\ U \\wl : ?{Q T ) = W U Wl sA (Q t ) + l|Vw|U s! (Q r ), 

H u II<^Wt) = W^wI^Qt) + \\V 2 u\\l s aQt) + liatu|k,,«M. 

Denote B(xo,r) the open ball in R 3 of radius r centered at xo and denote 
by B + (xo,r) the half-ball {x G B(xo,r) \ x$ > 0}. For zo = (xo,to) denote 
Q(z , r) = B(x , r) x (t - r 2 ,t ), Q + (z , r) = B+(x , r) x (t - r 2 , to). In this 
paper we shall use the abbreviations: B(r) — B(0,r), B + (r) = £> + (0,r) etc, 
B = B(1), B+ = B+(1) etc. 



2 Main Results 

Our work deals with the sufficient conditions of local regularity of suitable weak 
solutions to the MHD system near the plane part of the boundary. In [9] the 
following results were obtained. 

Theorem 2.1. There exists an absolute constant e* > with the following 
property. Assume (v,H,p) is a boundary suitable weak solution in Qt and 
assume zq — (xo,to) G dfl x (0, T) is such that xo belongs to the plane part of 
dft. If there exists ro > such that Q + (zq, ro) C Qt and 

^ J ( \v\ 3 +\h\ 3 +\ P \i y X dt < £*, 

°Q+(z .r ) 

then the functions v and H are Holder continuous on Q + (zo, ^). 

Theorem 2.2. For any K > there exists £o(K) > with the following prop- 
erty. Assume (v,H,p) is a boundary suitable weak solution in Qt and assume 
z o = (xo,to) G <9S1 x (0,T) is such that xq belongs to the plane part of <90. If 

limsupf - / |ViJ| 2 dxdt ) 1/2 < K (5) 

Q(z ,r) 
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/If \ V 2 

limsup - / |Vv| 2 dxdt < e , (6) 

r^o v r J ) 

Q(z ,r) 

then there exists p* > such that the functions v and H are Holder continuous 
on the closure of Q + (z , /?*). 

Let us explain the main differences between these theorems. The statement 
of the theorem 2.1 contains smallness conditions on the three functionals, but 
these conditions have to hold only for one value of cylinder radius. In theorem 
2.2 we have conditions for all sufficiently small values of radius, but smallness 
condition (6) is imposed only on the velocity v. 

To describe more conditions of local regularity we will need the following 
notations 

/ \ 1/2 

E(r) = ( i / \Vv\ 2 dxdt ) , 

V Q+(r) ' 

E*{r) = (l J \VH\ 2 dxdtV 2 , 

V Q+(r) 1 

A(r) = U sup / \vfdy) 1/2 , 

V te(-r 2 ,0)B+(r) 



A,(r) = U sup / \H\ 2 dy 

V te(-r 2 ,0) B +(r) 

C q {r)^(^ J \v\"dydt) 

V Q+(r) ' 

F q (r) = J iH^dxdt" 1 " 

V Q+(r) 

D(r)=(^ J \p - [p} B+ (r)\ 3/2 dydt) 



1/2 

1 

i/ q (7) 



2/3 



+ (r) 

, s 1.3 2/3 



D s (r)=Ri~U J ( / \Vp\ s dyY'~ 2 dt) 

— r 2 B+(r) 

C(r) = C 3 (r), F(r) = F 3 (r), D.(r) = Djg(r). 

Note that the equations (I), (2), as well as the functionals (7) and the state- 
ments of the previous theorems are invariant under the scaling transformations 



v P (y, s ) = p v (pv + xo, p 2 s + t ), 

H P (y,s) = pH(py + x Q ,p 2 s +t ), 
P P (y, s) = p 2 p(py + x , p 2 s + t ). > 



(8) 



We use the approach which was originally developed in [4] for the Navier- 
Stokcs equations (and later it was used also in [2]). According to this approach 
the regularity of solutions follows if one of the functionals (7) is bounded uni- 
formly with respect to r and additionally one of these functionals is small only 
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for a single sufficiently small value of the radius. Our goal is to obtain the same 
result for the solutions to the MHD system. 

The main result of our work is the following theorem, that is a kind of 
"interpolation" of theorems 2.1 and 2.2. 

Theorem 2.3. For arbitrary K > there is a constant £\(K) > with the 
following property. Assumw (v,H,p) is a suitable weak solution to the MHD 
system in Qt and z = (x ,t ) e dfl x (0,T) where x belongs to the plane part 

of an. if 

f 1 f \ 1 / 3 / 1 f \ 1 / 2 

limsupf — / M 3 dxdt j +{— / \H\ 2 dxdt J < K (9) 

Q+(za,r) Q+(z ,r) 

and one of the following conditions holds 

/ 1 [ \ 1 / 2 
liminfl - / \\7v\ 2 dxdt ) < Si, 
r^O \ r J ) 

Q+(zo,r) 

( 1 f \ 1 / 2 

liminf ( - sup / \v\ 2 dxdt ) < £\, 
v r - r 2 <t<0 J ) 



-r 2 <t<0 

B+(x ,r) 



( 1 f x 1 / 3 

liminf^ — / |w| 3 dxdt J < E\ 



(10) 

1/3 



Q+(za,r) 

then there exists > such that the functions v and H are Holder continuous 
on the closure of Q + (z 0} p* )- 

Note that it is possible to prove a lot of analogues of theorem 2.3. Generally 
the proof consist of two steps. The first step is the proof of boundcdness of 
energy functionals (7). Usually to do this it is sufficient to have boundedness 
condition for one functional depending on v and for another one depending on H. 
In our work this step is carried out in Section 4. Also we will use some estimates 
for the magnetic field H, that can be obtained if we consider equation (2) as 
the heat equation with lower order terms depending on v. This inequalities are 
proved in Section 3. 

The second step is the proof of regularity condition when all of functionals 
(7) are bounded and one of functionals on v is small for a single sufficiently 
small value of r. This result can be found in Section 5. 



3 Estimates of Solutions to the Heat Equation 

In this section we study solutions of the heat equations with the lower order 
tcrins' 

d t H - AH = div(v®H-H®v) in Q+ . 

U|x 3 =0 = 0, 

H 3 \ X3=0 ~0, H a ^\ X3=0 — 0, a = 1,2. 
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Namely, we assume the functions (v, H) possess the following properties: 

v, HeW£°{Q + ), 
^|x3=o = 0, ^31x3=0=0 in the sense of traces, 

for any r\ 6 Cq°(Q;M. 3 ) such that i]^\ X3=0 — the following integral identity 
holds 

J ( - H ■ d t n + VH : Vr^j dxdt = - J G :Vrj dxdt, ^ 
Q+ Q+ 
here G — v ® H — H ® w, and 

divw = 0, divii = a.e. in Q+. (13) 

Lemma 3.1. Assume that conditions (11) — (13) hold. Then for any < r < 1 
and < < 1 i/ie following estimate holds 

F 2 (9r) < c6 a F 2 {r)+cO-'iC{r)A*{r). (14) 



Proof. Denote by v* and H * the extensions of functions v and if from Q + onto 
Q. Fix arbitrary r <E (0, 1) and let ( G C°°(Q) be a cut off function such that 
C = 1 on Q(r) and supp£ Cfix (—1,0]. Denote II = R 3 x (—1,0) and denote 
by G the function which coincides with G* on and additionally possesses 

the following properties: G <G W-, ' (II) n L is e (II), G is compactly supported in 

11)5 

II , and 

l|G|U 62 (n) < c||G*|| L6aW(5)) < c\\G\\ Le 2(Q+(S)) (15) 

"5" ' £ ' "5" ' 

We decompose H* as 

if* = H + H, 



where if is a solution of the Cauchy problem for the heat equation 

( d t H — AH = divG in n, . . 

I H|t=-i=0, 1 j 

defined by the formula if = T * div G = — Vr * G, where T is the fundamental 
solution of the heat operator. The function if satisfies the homogeneous heat 
equation 

d t H-AH = in Q(§). (17) 
Take arbitrary 6 e (0, |). We estimate ||ff||z, 2 (Q+(0 r )) m the following way 

II^IU 2 (Q+(er)) < \\H*\\L 2 (Q(9r)) < \\H\\ L 2 (Q(0r)) + \\H\\ L2 {Q(6r)), ( 18 ) 

For \\ H \\L 2 (Q(er)) we have 

\\H\\L 2{Q{ er)) < c \\H\\ L2(Q{ r )y (19) 
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As H satisfies (17) by local estimate of the maximum of H via its L 2 -norm we 
obtain 



\\ H \\l 2 {Q(6t)) < C02 ||-ff|U 2 (Q(r)) < 
< C0i {\\H*\\L 2( Q ( r)) + \\H\\L 2(Qm ) 



(20) 



So, we need to estimate \\H ||l 2 (Q(§))- As singular integrals are bounded on 
the anisotropic Lesbegue space L s j (see, for example, [8]) for the convolution 
h = T * G we obtain the estimate 



\\M\wr (Q(r)) < c||G||i (n) . 

5' 2 5 ' 

On the other hand, from the 3D- parabolic imbedding theorem (see [1]) 

Wlt(Q)^W^(Q), as l-(| + f -|-|) >0, 
for p = q = 2 and s = |, I = 2 and for H = —Vh we obtain 



\H\\L 2 (Q(r)) < C \\G\\ Le my 

5'^ 



(Note that the constant c in this inequality does not depend on r). Taking into 
account (15) we arrive at 



\\H\\L 2 (Q {r)) < c ||G|| L| 
From the definition of G we obtain 



-(§))• 



(21) 



l G lli6 2 (Q+(§)) < c 

5 



\v®H\\l &(B+(r/2)) dt 



\-r 2 /4 



< 




V \\lB+(r)\\ H \\2,B+(r) dt 



< 



(22) 



< \\H\ 



2,oo,Q+(r) 



u \ ^ 

/ \Ht,B H r)dA <ric(r)A«(r). 



Combining inequalities (20)-(22) we will get the statement of lemma. 

□ 

Using interpolation inequality (24) for C(r) in the right hand side (14), we 
will obtain inequality (14) in another form 

Corollary 3.1. Assume that conditions (11) — (13) hold. Then for any < 
r < 1 and < 9 < 1 the following estimate holds 



F 2 (6r) < c9 a F 2 (r) + c0-? (r)Ai (r)A*{r). 



(23) 
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4 Boundedness of energy functionals 



In this section we derive estimates of energy functionals which allow us to ob- 
tain uniform boundedness (with respect to the radius) of all functionals (7) if 
boundedness of some of them is known. 

Observe that one can prove a group of estimates that are the consequences 
of Holder inequality, embedding theorem and interpolation inequality. 

C(r) < A^(r)Ei(r), F{r) < A} (r)[E} (r) + F 2 * (r)] (24) 

D(r) < c£>i(r), Di(r) < cD s {r), Vs > 1. (25) 
First of all we will prove the decay estimate for the pressure 

Lemma 4.1. Ifv,p,H are the suitable weak solution near the boundary to the 
MHD equations in Q + . Then for any < r < 1 < 9 < 1 the following 
estimate holds 

D12 (8r) < c9 a (Di2 (r) + E(r)) + 

, " i x (26) 

+c{6) [E{r)A^{r)C^{r)+E,{r)A}{r)F^{r)j . 



Proof. To obtain (26) we apply the method developed in [3], [5], see also [6]. 
Denote IL- = x (-r 2 , 0). We hx r £ (0, 1] and 9 G (0, ±) and define a function 
g : IT+ -> M 3 by the formula 

J rotH x H- (v ■ W)v, in Q+(r), 
9 - \ 0, in II+\Q+(r) 

Then we decompose v and p as 

V = v + v, p = p + p, 

where (v,p) is a solution of the Stokes initial boundary value problem in a 
half-space 

f d t v-M + Vp = g, + 
{ divw = r ' 

v\ t= o = 0, v\x 3 =o = 0, 
and (v,p) is a solution of the homogeneous Stokes system in Q + (r): 

U|x 3 =0 = 0. 

For Vj3 and Vp the following estimates hold (see [5], see also [7]): 



< 



ll V PlU 12 3 (Q+(r)) + 7ll V «IL 12 3 (Q+(r)) < 
TT ' 2 r TT ' 2 

c ( \\HxrotH\\ L 3 ( Q+(r)) + ||(«- V)«|| iia a (Q+(r)) ), 

\ 11 ' 2 11 ' 2 / 



(27) 
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l|Vp|| iM ,(0+(flr)) < Cfl° ( -||V«|| Ll , 3(Q+(r)) + ||Vp|| Lia 3 (Q+(r)) )• (28) 

11 ' 2 V / 11 ' 2 11 ' 2 / 

To estimate the right hand side of (27) we will use Holder and interpolation 
inequalities 




Vv\\l\\v\\idt\ 



< ||V«|| 2 I / IMM I < ||Vd| 2 I / \\v\\l\\v\\Ut I < 



(29) 



< I|VW||2||V||| I00 ||«||| 

Term \\H x rotH\\ Ll2 3 (Q+(r)) can be estimated similarly. 

TT'2 

The right hand side of (28) can be estimated as follows 

-||Vw||l 12 3 (Q+(r)) + l|Vp|| Ll2 3 (Q+(r)) < 
r TT' 3 TT'3 

< c(||Vu|| 2 ,Q+( r ) + \\Vp\\ Ll2 s(Q+(r))+ (30) 

\ TT ' 2 

+ 7ll Vfi llL 12 s(Q+(r)) + l|Vp|| Ll2 3 (Q+(r)) 
r TT ' 2 TT ' 2 

Combining inequalities (27)-(30) we will get the statement of lemma. 

□ 

Theorem 4.1. Let v,p,H are the suitable weak solution near the boundary to 
the MHD equations in Q + and 

C{R) + F 2 {R) < M, < R < 1. 

If we consider the following functional 

25 

£ (r) =A 2 (r)+E 2 (r) + A 2 (r) + E 2 (r) + £>g (r) , 

ii 

then the following estimate will hold 

L(r) <C(M)(r a L(l) + l). 

Proof. From local energy inequality we obtain 

L{9r) < c (Cl{26r) + F 2 2 (26>r) + C 3 (28r) + C(26r)D{29r) + 

C 2 (26r)F 3 (28r) + C(28r)F 2 (29r) + D]l(8r)) . 

ii / 



Now we will estimate every term in the right hand side. Our goal is to prove 
the following estimate 

£{6r) <-L(r)+C(M). (31) 

Then we can use a standard iteration procedure (see [4]) and obtain the state- 
ment of the theorem. 

Estimates for the first three terms are obvious. To estimate the 4th term we 
use Young inequality and inequality (25) 



C(26r)D{29r) < c (£>§(26>r) + M 25 ) . 



Now we are going to prove estimate for Dr2.{29r). To do this, we will use 
inequalities (24) and (26) 

D^(26r) < 

< c6 a (D 12 (r) + £7(r)) + c(0) (E(r)A* (r)C^ (r) + E* (r)A? (r)F 3 * (r)) < 
< c9 a (D^(r) +E(rfj + c{6) (L*{r)M% + L% (r)Kf (rfj . 

(32) 

We use interpolation inequality to estimate F 3 (r) 

F 3 (r) < Ff (r)Fj(r) < c (aI(t) f^|(r) + F$ (r))) 6 F« (r) < 

- - v v yy (33) 

Now we substitute this to (32) 

£>i2(20r) < c6» Q (r>i|(r) + E(r^j + 
+c{6) (£,i(r)Mi +L^(r)M^ +L§( r )M^ . 

As the result we obtain 

£>jf (2<9r) < c0 Q £(r) + c(6) [L^{r)M kl + £ts (r)M fea ) . 

Since the right hand side of the last inequality contain £ (r) in the degree smaller 
then 1, choosing 9 sufficiently small and using Young inequality we obtain an 
estimate (31). 

To estimate the last two terms we use (33) 



C(r)F 3 2 (r) <c(zi{r)M^ + Z^Mi 



and Young inequality. The second term can be estimated in the same maner. 

As the result we obtain (31). Next by standard iteration procedure we finish 
the prove of the theorem. 

□ 
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5 Proof of main results 

As a first step we obtain theorem 2.1 without smallness condition on a pressure. 

Lemma 5.1. For arbitrary M > there is £\(M) > 0, such that ifv,p,H are 
the suitable weak solution near the boundary to the MHD equations in Q + , 

A(R) + E(R) + A*(R) + E*(R) + F 3 (R) + D^(R) < M, V0<R<1 (34) 

and 

C(l)+F 3 (l)<ei, (35) 

then the functions v and H are Holder continuous on Q + (r tt ) for some < r* < 
1. 

Proof. Assume that the statement of the lemma is false. Then there are 
sequences of v n ,p n ,H n of suitable weak solutions in Q + , such that 

C(v n , 1) + F 3 (H n , 1) = e n -> 0, asn^oo (36) 

and is a singular point. Then by theorem 2.1 

C{v n ,r) + D{ Pn ,r) + F 3 {H n ,r) > e* (37) 

for all < r < 1. 

On the other hand from (26), (34), (36) and the embedding theorem we have 
D(p n ,r) < cD^{p n ,r) < cr a M + c{r)M^e\. (38) 
So we fix < r < 1 and pass to the limit by n in (37) and (38) 
£* < lim sup (C(v n ,r) + D(p n ,r) + F 3 (H n ,r)) = 
= lim sup D(p n ,r) < cr a M. 

n— >oo 

As the result we obtain, that the inequality 

e* < cr a M 

must be true for arbitrary < r < 1. So we have a contradiction. 

□ 

Theorem 5.1. For arbitrary M > there is £2{M) > 0, such that ifv,p,H 
are the suitable weak solution near the boundary to the MHD equations in Q + , 
satisfying to (34) and one of the following conditions holds 

E(l)<e 2 , (39) 
A(l) < e 2 , (40) 
C(l) < £ 2 , (41) 

then the functions v and H are Holder continuous on Q + (r*) for some < r* < 
1. 
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Proof. The proof of this theorem is similar to the proof of previous lemma. 
We begin from the case (39). Let v n ,p n ,H n are the sequences of suitable weak 
solutions to the MHD system, such that (34) holds, 

E(v n , 1) = e n -> 

as n — » oo, and z = is a singular point. Then from lemma 5.1 we have 

C(v n ,r)+F 3 {v n ,r)>e 1 (42) 

for arbitrary < r < 1. 
On the other hand 

C(v n ,r) < 4-C(«„,l) < 4-^(« n ,l)^'(«n,l)-^0 (43) 

7" 3 ra 

as n — > oo and for any fixed < r < 1. From (23) we obtain 

lim sup F 2 (#n, r) < cr a M. (44) 
Next we use interpolation inequality 

^3 (H n , r) < F* (£T n , r)F J (H„ , r) . (45) 

3 

To estimate the second factor in the right hand side of (45) we use (24). So 
from (42)-(45) we obtain 

£i < lim sup (C(w„,r) +F 3 (u„,r)) < cr Ql M fc VO < r < -, 

n— >oo ^ 

and, if we choose r sufficiently small, we will have a contradiction. 

Observe, that E(r) and A(r) take part in (23) symmetrically, so the proof 
of this theorem in the case (40) is similar to the previous one. In the case of 
(41) for obtaining (44) is sufficient to use (14). 

□ 
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